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Abstract: The reaction process A + B — > is modelled for ballistic reactants on an 
infinite line with particle velocities va = c and Vb = — c and initially segregated conditions, 
i.e. all A particles to the left and all B particles to the right of the origin. Previous 
models of ballistic annihilation have particles that always react on contact, i.e. pair-reaction 
probability p = 1. The evolutions of such systems are wholly determined by the initial 
distribution of particles and therefore do not have a stochastic dynamics. However, in this 
paper the generalisation is made to p < 1, allowing particles to pass through each other 
without necessarily reacting. In this way, the A and B particle domains overlap to form a 
fluctuating, finite-sized reaction zone where the product is created. Fluctuations are also 
included in the currents of A and B particles entering the overlap region, thereby inducing 
a stochastic motion of the reaction zone as a whole. These two types of fluctuations, in the 
reactions and particle currents, are characterised by the intrinsic reaction rate, seen in a 
single system, and the extrinsic reaction rate, seen in an average over many systems. The 
intrinsic and extrinsic behaviours are examined and compared to the case of isotropically 
diffusing reactants. 
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I. INTRODUCTION 



A detailed understanding of reaction systems is an essential ingredient for the study of a 
broad range of problems M. In systems of many interacting reactants, it often happens 
that knowledge of the precise physical mechanism whereby reactions occur is irrelevant in 
determining the macroscopic behaviour. The most important factors are the number of 
different reacting species that combine in a single reaction, the type of motion that each 
of these species performs and the initial/boundary conditions in the vessel containing the 
reactants. Thus, simple theoretical models that mimic these elements can be used to study 
real systems as diverse as exciton dynamics in polymer chains [§] , monopole annihilation in 
the early universe || as well as the more conventional chemical processes M . However, even 
such highly simplified models can prove difficult to analyse. This is especially true if the 
reactions are not fully reversible as the methods of equilibrium statistical mechanics cannot 
be used, due to the lack of detailed balance. 

Many different, idealised reaction systems have been studied analytically over the last 
few decades Jj] - |22| usually for the case of fully irreversible reactions. These include same 
species (nA — > 0) and multi-species (J2k n kAk — > 0) processes, where n k of each distinct A k 
species combine to form a single inert product 0. The multi-species reaction systems exhibit 
particularly rich behaviour because reactions can only occur at places where all the necessary 
constituents are present - so called reaction zones. Often, due to the initial/boundary 
conditions or spontaneous symmetry breaking ([I3|,|]14|) each species of reactant is largely 



confined to its own domain. The places where these domains overlap usually takes up only 
a small part of the whole system, often causing the net reaction rate to differ drastically 
from that predicted by the mean-field-like rate equation. Such behaviour can give rise to 
complicated structures, and is of particular interest in the context of pattern formation and 
growth determination in organisms |[23|| . The physics of multi-species reaction processes is, 
therefore, largely determined by activity in and around the reaction zones. The important 
factors being the way reactants flow from their respective domains into the overlap region 
as well as the behaviour of the reaction mix inside the reaction zone. 

In this paper a model of the two-species system A + B — > is introduced for the case 
of ballistically moving reactants and a pair-reaction probability less than one. The model 
is then solved exactly, allowing an analytic study of the dynamic reaction zone formed at 
the overlap of the A and B domains. However, before describing the model in detail two 
important statistical quantities associated with the reaction zone are defined. 
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The intrinsic and extrinsic reaction rates 



A comprehensive analysis in the RG framework for the two-species A + B — > reaction 
system with diffusive reactants [O identified two distinct sources of fluctuations affecting 



the behaviour of the reaction zone, first inside the reaction zone itself and second in the 
currents of particles entering the reaction zone. 

The shape and size of the reaction zone is determined by the typical lifetimes and motions 
of reactants in the overlap region. These lifetimes in turn depend on the fluctuating local 
densities of reactants, as well as the reaction probability for an interacting pair of A and B 
particles. As the density fluctuations will be caused (in part) by reactions that have already 
occurred, the reaction rate can become highly correlated both spatially and temporally. This 
behaviour in the overlap region leads to the following definition for the intrinsic reaction 
rate. 



• The intrinsic reaction rate lZi(x r ,t) measured at a distance x r from the centre of the 
overlap region at time t, is the typical spatial reaction rate, i.e. production rate of 
particles per unit time per unit space, seen in a single realisation of the system's 
evolution. The intrinsic reaction rate characterises the intrinsic reaction zone - the 
instantaneous reaction zone formed where the A and B domains overlap. 



A second source of fluctuations can come from the A and B currents that flow into this 
(intrinsic) reaction zone. Even if these currents are (on average) equal in magnitude, any 
fluctuations about the average will contribute a noisy component to the reaction zone's 
motion, i.e. it will move stochastically about its expected position. Hence, given some 
initial conditions only a probabilistic statement can be made about the reaction rate at 
some later place and time (x,t). This suggests the definition of a second quantity, the 
extrinsic reaction rate. 



• The extrinsic reaction rate TZE{x,t) is defined as the probability density for reactions 
to occur at a time t and position x given some initial distribution function for the 
particles' positions. Therefore, TZe is the expected reaction rate found by averaging 
over all allowed realisations of the system's evolution. 



Models with ballistically moving reactants have been extensively studied in the context 
of the reaction kinetics in an ideal gas, where the mean free path of the reacting particles is 
similar to their separation 0-0. In particular the two-species case, with particles A and 
B having velocities v a = +c and vb = — c respectively, has been studied with homogeneous 
initial conditions, i.e. each species initially randomly placed throughout the infinite line. 
More recently, the initial conditions of segregated particles was studied (P],[^Q[) with the As 
initially to the left and the Bs to the right of the origin, again with reactants that always 
annihilate on contact. The initial positions for each species of particle (confined to their 
own domains on either side of the origin) were chosen to be random, thereby introducing 
fluctuations into the particle currents flowing into the reaction zone. By averaging over 
all initial distributions the form of the extrinsic reaction rate was derived and shown to be 
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Gaussian. However, because there can be no reaction fluctuations in this model (the particles 
always annihilate on contact) the A and B domains never overlap and the instantaneous 
reaction zone has effective width zero. Also, as the reactants move ballistically the system 
does not have a genuinely stochastic evolution, i.e. once the initial conditions (the particle 
positions) are fixed so is the system's future. 

This model has a generalisation to the case of arbitrary reaction probability; this is the 
general model introduced below. The fluctuations in the particle streams can be retained 
and combined with the fluctuations in the (now) finite-sized overlap region. This then 
allows the intrinsic and extrinsic behaviour of the A + B — > reaction system to be studied 
analytically. 

The definition of the model 

The model consists of a one-dimensional continuous space in which two species of reactants, 
A and B particles, move with fixed velocities v\ = +c and vb — —c, see figure 1. Initially, 
the reactants are separated, i.e. at time t — A particles are distributed in the interval 
[— oo : 0] at positions (yi,U2 • • •) and the B particles are in the interval [0 : oo] at positions 
(zi, z 2 ■ ■ •)• The subscripts on y m and z n refer to the relative initial order of the particles 
counted from the origin. Because of the ballistic motion of the reactants, the trajectories of 
the particles retain their initial ordering for all time. 

Two different distributions for {y} and {z} will be considered, first equally-spaced re- 
actants (section III) and second, random positions of the reactants (section IV). For both 
cases the average density is chosen to be g, leading to average particle currents of ±cg. 
However, the initial conditions studied in section IV will be shown to introduce Gaussian 
fluctuations about these average values. 

The position x and time t that the mth A and nth B particles' trajectories intersect are 

x = (y m + z n )/2 t = (z n - y m )/2c (1) 

When such an (m,n) pair of reactants' trajectories meet, there are three distinct events 
that can occur. 

• If both reactants still occupy their trajectories a reaction can occur with probability 
p. If a reaction occurs both particles are removed from the system and a product 
particle is considered to have been deposited at the point of annihilation. The merely 
serves as a marker and plays no further role in the evolution of the system. 

• If both reactants are still travelling along their trajectories then with probability q = 
(1 — p) no reaction occurs and the particles continue unaffected. 

• If one of the particles has previously been annihilated, i.e. one trajectory is unoccu- 
pied, the other particle continues unaffected with probability 1. (Of course, if neither 
particle is present no change occurs when the empty trajectories cross.) 
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Thus, two distinct sources of fluctuations are included in this model, in the reactions (if 
p < 1) and in the currents of particles flowing into the reaction zone (if {y} and {z} have 
random elements). Hence, the system models a finite-sized stochastically moving reaction 
zone, with measurable intrinsic and extrinsic behaviour. However, it is clear that this model 
represents a special case, in that the two sources of noise are uncoupled. The fluctuations in 
the trajectories are quenched at t = 0, and the probability that an (m,n) pair annihilates 
depends only on the total number of trajectories crossed by each particle, and not the 
trajectories' exact positions. This uncoupling of fluctuations allows the probability density 
V for an (m, n) pair to react at (x, t) to be written as the product of two independent 
distributions 

V(m,n,x,t) = P(m,n)G mn (x,t) 

where P(m,n) is the probability for an (m, n) pair to mutually annihilate and G mn {x,t) 
is the probability density for the trajectory intersections ([!]). The calculation of the total 
annihilation probability can therefore be decomposed into, (i) a counting problem for the in- 
teger variables (m, n), and (ii) the derivation of the distribution functions for the continuous 
random numbers {y} and {z}. 

The remainder of the paper is structured as follows. In section II the probability of 
pair annihilation P(m,n) is derived ([12]) by mapping the model onto a simple system of 
target particles on a one-dimensional lattice. The behaviour of the lattice system is then 
briefly examined in the context of radiation damage of crystals. In section III the form 
of the (intrinsic) reaction rate for the simple case of equally-spaced reactants is studied. 
In particular, the steady-state reaction-rate (|T6|) and particle densities (|17D, and the time- 



dependence of the reaction rate fllTf ) are derived. Finally, in section IV fluctuations in the 
initial particle positions are treated. The distribution function Q mn (x, t) is found and used 
to derive the forms of the intrinsic and extrinsic reaction rates, (|24]) and (|25|) respectively. 
The appendix shows how the calculation in section II may be translated into the second- 
quantisation formalism and relates the mapped system to the algebra SU q {2). 



II. THE PAIR-REACTION PROBABILITY 

The main result of this section is the derivation of the (m, n) pair annihilation probability 
P(m,n), equation (|12"D. In section I, it was noted that this quantity is independent of the 
initial positions of the particles, depending only on their relative order. Making use of 
this fact, a simpler system can be treated, that still preserves the order that the A and B 
particles pass through each other. 

Consider now a one- dimensional, semi- infinite lattice (with sites n = 1,2---) with a 
B particle initially occupying each site. At discrete time steps (m = 1, 2 • • •) the mth A 
particle is 'shot' through the B array, passing sequentially through each site until its eventual 
annihilation, see figure 1. An A particle either moves through a site with a B present with 
probability q, or reacts with a B at that site with probability p. Once an A and B pair have 
reacted, the site that the B occupied becomes vacant and any subsequent A that passes 
through the vacant site does so with probability 1. This whole process, of a single A passing 
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through the lattice and eventually reacting, is considered to happen in a negligible amount 
of time. This mapped system is merely a deformation of the coordinates used in the original 
model, so both systems share the same pair annihilation probability P(m,n). 

It is useful to consider the statistics of the positions of the m vacancies that exist in the 
lattice immediately after the mth time step. The positions of these vacancies are labelled 
(ni---n m ), where it is important to note that the subscript used refers to the relative 
positions of the sites, i.e. n\ < n 2 < ■ ■ ■ < n m , and not which A particle caused the vacancy 
at that position. 

The rest of this section is devoted to the calculation of P(m, n), the probability that the 
mth A reacts at site n. First, the probability of a particular distribution \l/ m (n 1 ■ • -n m ) for 
the positions of the m vacancies is derived, equation (|3|). These probabilities can be used 
as a convenient basis, in the sense that all other probabilistic quantities may be expressed 
as linear combinations of the {\F m }. This basis is then used to find the expected vacancy 
density V m {n) at site n, after m A particles have passed, equation (|9]). As an aside, a parallel 
is drawn between this simplified model and a crystal that has been damaged by radiation. In 
particular, it is shown that the damaged region described by V m {n) propagates like a soliton 
through the B array, equation fllPl) . Finally, the discrete gradient of V m {n) is then used 
to calculate the required quantity P(m,n), equation ([T2|) . The method described below 
translates into the second-quantisation formalism and shows the system to be described 
by the algebra SU q (2). The techniques used in this formalism are briefly reviewed in the 
appendix. 

The basis for the mth time step ^> m 

The probability \l/ m (ni • • -rz m ) that after m A particles have passed through the B lattice 
vacancies exist at sites n\ <n 2 < ■ ■ ■< n m is now derived. Consider first the simplest case 
m=l, i.e. just after the first time step. The probability that the A particle has annihilated 
with a B on site n\ producing a vacancy there is 

tti(ni) = V q ni ~ l = (q' 1 - i)q ni 

where p is the probability a single reaction could occur, and q = (1 — p). 

Now consider the state of the system after the second time step, i.e. after a total of two 
A particles have passed through the B lattice. Vacancies now exist at sites ni and n 2 (where 
relabelling may be necessary to ensure that n\ < n 2 ). There are two histories that contribute 
to this configuration, each with different probabilities. Either a vacancy first appeared at 
site n 2 and then the second one at site ni, or the first vacancy at site ni and the second at 
site n 2 . The probability of the first history is just the product ^i(ni)^i(n 2 ) as the second 
A particle does not pass through the vacancy produced by the first A particle. However, in 
the second history described, the second A particle does pass through the vacancy at 
increasing its chance of reacting with any site n > n\ by a factor of q~ x . Therefore, 

* 2 (n 1; n 2 ) = *i(rii)#i(n 2 ) + g _1 *i(ni)*i(n 2 ) 

# 2 (ni,n 2 ) = (1 + <T 1 ) 1 I'i(rii)^i(ri 2 ) (2) 

Mnnn 2 ) = (q~ 2 - l)^" 1 - l)q n ^ . 
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It is simple to generalise to the m- vacancy basis ^ m {n\ ■ ■ ■ n m ). The states at the (m — l)th 
time step that can contribute to an mth time step configuration are those with vacancies 
at (n 2 ---n m ), n 3 ■ ■ ■ n m ), (ni ■ ■ ■ n m _i). In each of these cases, to produce the 
final state (ni---n m ) the mth A particle must react with sites ri\,n 2 - ■ -n m respectively. 
Therefore, taking account of how many vacancies the mth A particle must pass through in 
each case, a relation between the bases that describe the system after time steps (m — 1) 
and m can be written 



^ mijll ' ' ' rim) 



' rim ) 

+q~ 1 ^ 1 (n 2 )^m-i(rii, n 3 • • • n m ) + 
+q- ( - m - 1 ^ 1 (n m )^r 



_i(ni,n 2 • ■■n m -i / 



Assuming that it is possible to write ^ m -i = C m ~ 1 ITj=i ^i( n j) an d using the result (fj), 
the form for general m follows by induction 



[Tli 



^m(ni 



n r . 



rir, 



l + g- 1 + ---g- (m - 1) )C m _ i n*K) 



C m II *(«j0 
3=1 



n 

j=l L 



q-' - 1 
q- 1 - 1 



n[(?- 



(3) 



3=1 



Hence, the distribution of the m vacancies is given by a product of exponentials in the site 
labels, with care being taken to preserve the order n\<n 2 <- ■ ■< n m in any sums that they 
appear. 



The vacancy density V m (n) 

The probability that a vacancy is found at site n after m A particles have passed through 
the lattice, V m (n), is now derived. This density is given by the sum of all the \l/ m that 
include a vacancy at the site n, 

m 

k=i 

where the notation \1/ is used to denote an internal sum over all the unfixed variables, i.e. 
any rij with j ^ k, that respects the order rii<n 2 <- ■ -<n m . Hence, 

§ m (ni • • • ,n k = n, ■ ■■n m ) = 

n—(k— 1) n—1 oo oo 

"•• J2 J2 ■■■ J2 y m (ni---,n k = n,---n m ). (4) 

"1=1 n fc _!=n fc _ 2 +l n k+1 =n+l n m =n m _i+l 

For the case m — 1, the vacancy density is simply V\(n) = ^ r i(n). However, after the second 
A particle passes through the B array a sum must be made over the unfixed variables 

V 2 (n) = '^ 2 (n 1 = n,n 2 ) + ^ 2 (n 1 ,n 2 = n) 

oo n—1 

v 2(n) = *2(ni = n,n 2 ) + 51 *2(ni,n 2 = n) 

n2=n+l ni=l 
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( 00 71— 1 

j=n+l j=l 

(oo 
£>i(;)-*i(n) 

V 2 (n) = g n (<r 2 -l)(l-^i(n)). 

The density for m = 2 is therefore related to the m = 1 case, by making use of the product 
form of ^2- It is possible to generalise this result and obtain a recursion relation. The 
following two results will prove useful. First the product form of \l/ m is used to relate v P m 

to ^ m _i 

ty m (ni---,n k = n,---n m ) = q n (q~ m - l)* m _i(ni • • •n k - 1 ,n k+1 ■ ■ -n m ). (5) 
Second, a slightly less trivial result, that nevertheless has a simple interpretation 

m m— 1 

= 1 - E * m _ 1 (n 1 • • -,n k = n, ■ • -n m _i) (6) 
fc=i 

= 1 - V m -An). 



*m-i(wi • • •n fc _i,7i fe+1 • • -n r , 



fc=i 



The LHS of this equation is the sum of all the ^ m _i that do not include a vacancy on site n. 
This is simply the sum of all possible ^ m -\ (= 1 by normalisation) less those that include 
the site n; the statement on the RHS of (||). The \Ps in this equation are identical to those in 
equation (|J) in as much as they involve an ordered sum over all unfixed variables. However, 
it should be noted that, though on the LHS the term n k = n has been factored out by using 
(|), the order restriction still holds, i.e. n k -i < n and n k+ i > n. 

Both these results generalise the method used already to find V^n). Therefore, for the 
case of m vacancies 



V m (n) = ^2^f m (ni---,n k = n,---n m ) 

k=l 

m 

V m {n) = q n (q- 1 
V m (n) = q n (q-' 



1) E *m-i(ni • • ■n k _ 1 ,n k+1 ■ ■•n m ) 

k=l 

l)(l-y w -i(n)). 



(7) 
(8) 



The recursion relation (|j) can be solved, with the initial condition Vi(n) = ^i(n), to give 
the vacancy density at site n after m A particles have passed through the B array 



V m (n) 
V m (n) 



E 

3=1 



-q 



k=j 

m— 1 
(n— to) ^ ' 

3=0 



n—k\ 



(9) 



(n-m)j 



fc=0 



Interestingly, in the limit of large m (in particular q m <C p) the vacancy profile depends 
purely on the difference s = (n — m). Interpreting the system as a crystal (the B array) 
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that is being damaged by incoming radiation (the A particles), this implies that the interface 
between the damaged and undamaged regions reaches a steady-state moving profile 



lim V m (n) = V(s) 



<z s £(-<z s )V (i+1)/2 

j=0 

-{1- tanh(ps/2)) + 0(p) 



(10) 



(11) 



where the second equation becomes valid in the case of low reaction probability p <^ 1 and 
is analytically continued for s < 0. Therefore, the damaged region propagates like a soliton 
through the B array. 



The pair-reaction probability P(m, n) 

The difference between the vacancy density at site n just after the mth and (m — l)th time 
steps is the probability that at the mth time step an annihilation occurs at site n. This is 
the required pair reaction probability P(m,n). 



P(m, n) 



V m (n) - V m -i{n) 



P(m,n) = q {n ~ m) £ 



m—l 



3=0 



k=l 



(12) 



which is the main and final result of this section. For the case studied in 0, i.e. q — 0, this 
result reduces to a delta function as expected, i.e. P(m,n) = 5(m,n). In the next section 
P(m, n) is used to study the original model described in section I, with equal spacing between 
neighbouring reactants at t — 0. 



III. THE REACTION ZONE WITHOUT CURRENT FLUCTUATIONS 

Reverting back to the original model described in the first section, P(m, n) is reinterpreted 
as the probability that an (m, n) pair of particles mutually annihilates in a system of bal- 
listically moving reactants. In this section, the model is studied with the initial conditions 
y m = —mg" 1 and z n = rig' 1 , i.e. with no fluctuations in the particle currents. The point of 
intersection of an (m, n) pair of trajectories is therefore 

x = x r = (n — m)/2g t = (m + n)/2cg, (13) 

i.e. reactions occur only at discrete positions and times. The lack of fluctuations in the 
currents negates the need to discuss the extrinsic behaviour, as the intrinsic reaction zone 
will not wander stochastically. In fact as TZi = TZe for these initial conditions, the zero 
current-fluctuation reaction zone will be denoted by IZq to avoid confusion with section IV. 
The steady currents also mean that the centre of the overlap region is always at the origin, 
and therefore x r = x in this case. As will be shown in the next section, even when current 
fluctuations are taken into account, the form of the intrinsic reaction rate remains the same 
as for this CclS6 ; clS lone; then measured from the centre of the stochastically moving 

overlap centre seen in section IV. 
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In the rest of this section several aspects of equation fll2|) are studied. First, the steady- 
state reaction rate flllf ) and particle densities are derived flT7p. Second, the dynamics of 
P(m,n) in the 0(p) approximation is studied and the system is shown to decay to the 
steady state in finite time (|T^). Finally, the correlations in the reaction zone are briefly 
discussed, with reference to the mean-field and 0(p) approximations. 



The steady-state limit and 0(p) time dependence 

After a sufficient length of time, characterised by q m -C p, equation ([12]) becomes a function 
of (n — m) = s only. This is in much the same way as for equation ([[(]) and can also be 
interpreted as a late-time or steady-state limit. Therefore, for large m (or n by symmetry) 
any (m + I, n + I) pair will have the same probability of reacting, regardless of the value / 
takes, 

oo 

lm^ P(m,n) = P{s) = q s £ [(-q s ) j q j(J+1)/2 (l - q 1+1 )] (14) 



3=0 

The production rate of particles at a position x r = s/(2g) can now be found, namely 
2cgP(2x r g). However, the nature of the initial conditions means that reactions can occur 
only at discrete positions (s — 0, ±1 • • •). So, to derive the reaction rate TZ (x r ) which is the 
production rate per unit space, a coarse-graining is made over a length scale g^ 1 

n (x r ) = 2cg 2 P{2x r g). (16) 

This reaction rate (characterising the intrinsic reaction zone) is shown graphically in figure 
2 ior p = q = 1/2 with comparison made between the O (p) and mean- field approximations 
(to be described below). The particle densities can also be found, given that [1 — V(s)] is 
the probability a B particle still exists just after a reaction occurs at site s (and similarly 
for the A particles). Therefore, the density profiles gA{x r ) and gB(x r ) in the steady-state 
are 

g A {x r ) = g[l -V(-2gx r )] and g B (x r ) = g\l -V{2gx r )}. (17) 

These profiles are also shown, for the same parameters q = p = 1/2, in figure 2, with 
comparison made to the 0(p) approximation. 

To examine the passage to the steady state, it is illuminating to derive the time depen- 
dence of (0) to first order in p, i.e. take the limit p 1 in equation (|T2"D but still keep m 
finite. 

P(m,n) = p(2cosh{p(n-m)/2) -exp- p{m+n)/2 y 2 + 0{p 2 ) (18) 
-o l +\ pcg 2 8{cgt - 1 - g\x r \) 2 

K (x r ,t) = — ——2 +0(p), (19) 

2 {cosh{pgx r ) — exp pcQt ) 

where 6{u) = 1 or for u > or u < 0, respectively. It is clear that the system reaches a 
steady state exponentially quickly (a characteristic of systems with short-range correlations) 
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with a decay time (peg) 1 . To leading order in t, this exponential decay time will also be a 
feature of the model with randomly distributed particles to be examined in section IV. 



Correlations, the 0(p) and mean-field approximations 

The first order in p approximations (|TT|) , ( |15|) , (|18D and (|l9| ) valid for p <C 1 represent 



the case where each particle passes through very many of its opposite kind before finally 
annihilating. In this p <C 1 limit it could therefore be expected that a mean-field approach 
becomes exact. For the reaction model the mean-field approximation is one that neglects 
density correlations in the reaction zone. It assumes that the probability of a reaction 
occurring depends on the product of the probabilities that an A and B particle are present, 
i.e. that the A and B densities are independent statistical quantities. Therefore, given that 
[1 — V m -\ (n)] is the probability that the nth B has not annihilated before meeting the mth A 
particle (and similarly for the mth A particle by symmetry), the mean-field approximation 
is 

P m f(m,n) = p[l - K-i(m)][l - V m -i{n)\ 

P mf (s) = p[l-V(s-l)][l-V(s + l)} (20) 

where the second equation is the steady-state limit. From these equations it can be seen 
by substitution that the 0(p) approximations (valid only in the limit p — > 0) are indeed 
mean-field like. Though, this is not to say that the 0(p) and mean-field approximations are 
equivalent outside the small p limit, as can been seen in figure 2 where comparison can be 
made for p = 0.5. 

For the p < 1 case, the reaction zone will be much larger than the interparticle spac- 
ing, and any structure will slowly vary with the variables m or n. Therefore, the 0(p) 
approximations can also be captured from the continuum mean-field equations 

9qa dg A dg B dg B 

-£T = -c-z tqaQb -^rr = rg A g B . 

at ox at ox 

where r is the reaction parameter. However, the full results (H), (|10D , fll2"|) and (|T4]) cannot 
be obtained by such a mean-field approach, as can be seen in figure 2 where P m f(s) and 
P(s) are compared. The figure shows that the correlations are strongest in the centre of the 
reaction zone and decay quickly towards the edge of the overlap region. 

Interestingly, a similar cross-over to mean-field behaviour is also seen in the case of diffu- 
sive reactants Jl2[ in the limit of high diffusivity and low reaction rate. The interpretation, 
of many interactions between particles before final annihilation occurs, is the same. 
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IV. REACTION ZONES 



In this section the case of fluctuations in the initial conditions is treated (though still with 
each species initially segregated to either side of the origin). This allows the intrinsic 
and extrinsic reaction rates for the system with fluctuations in both reactions and particle 
currents to be found, equations (^4j) and (|25|) . 



The initial positions of each species of particle are chosen to be uncorrelated with an 
average spacing g" 1 and are described by Poissonian distributions with the variances of the 
interparticle spacing being g~ 2 . Hence, both y m and z n , the initial positions of the mth A 
and nth B particles, can be written as sums of independent random numbers 

Vm = (Vm - Vm-l) + (Vm-1 ~ Vm~2) H V (V2 ~ Vl) + {Vl ~ 0) 

Z n = {z n — Zn-l) + {z n -l — Z n - 2 ) H h fa ~ Zi) + (zi ~ 0) . 

The central limit theorem can now be used, stating that for large to and n such sums of many 
independent random numbers, drawn from the same distributions, become the Gaussians 

y{y m ) oc exp I — I Z(z n ) oc exp I — I . (21) 

The probability density for the mth A and nth B trajectories to intersect at position x 
and time t (the function Q mn {x,t) defined in section I) is now easily calculated. As both 
x and 2ct are linear combinations of independent Gaussian variables, for s m their 
distribution functions are also Gaussians with means and variances linear combinations of 
those in equation (fill) 



Gmn(x,t) = ^ -exp 

7r{n + m) 



(2gx — (n — m)) 2 (2cgt — (m + n)) z 
2(m + n) 2(m + n) 



(22) 



Before calculating the intrinsic and extrinsic reaction rates for p < 1, the deterministic 
case p = 1 examined in is first reviewed. Because in this case, the reaction probability 
P(m, n) = 5 mn the two domains of A and B particles never overlap and only an extrinsic 
reaction rate can be meaningfully defined. The probability density for reactions to occur, 
i.e. for (m, to) pairs to intersect, at (x,t) is given by 

Kp=i - 

K P =i - 

in the limit of large m. This is the extrinsic reaction rate for the p — 1 model and it implies 
the reaction front is a Gaussian random walker, covering a typical distance ~ (ct/g) 1 ^ 2 in a 
time t. 

For the general case p < 1 the points of intersection of these (to, m) pairs will be used as 
a convenient definition for the centre of the stochastically moving reaction zone. The forms 
of the intrinsic and extrinsic reaction rates can now be calculated. 



(m,n) 



X 



m=l n=l 



eg 



^/TTCgt 



exp 



{gxf 
cgt 



(23) 
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The intrinsic reaction rate 



The intrinsic reaction rate as defined in section I, is the reaction-rate profile seen in a single 
realisation of the system's evolution (measured relative to the centre of the A and B domains' 
overlap region). In this section it will be argued that, if reactions are measured relative to 
the position of the (m, m) pair intersections, the intrinsic reaction rate is on average equal 
to TZq, the reaction rate for the zero current-fluctuations case examined in section III. 

Consider the set of intersecting particle pairs that can be written (m — k,m + k) with m 
fixed and k varying. These pairs all share the same average time of intersection. Therefore, 
it is convenient to define the relative coordinates of intersection (x r , t r ) of these pairs to the 
central (m, m) pair 

2x r {jjm—k Zm+k) (l/m ^m) 2ct r ^m+k Vm—k) Vm) 

Xr = kg- 1 + Oik^g- 1 ) t r = + 0(k l/2 (2cg)- 1 ) 



where the 0(k l l 2 ) deviations are those expected from the Gaussian fluctuations ([H]). As 
the deviations are not functions of m, the noise in the particles' initial positions does not 
produce any time-dependent dispersion, i.e. the statistics of the pair intersections in the 
overlap region reach a steady-state. The relative time and positions are on average the same 
for the equally-spaced case seen in section III, though the Gaussian fluctuations about these 
average values introduce a ~ k x l 2 g~ x uncertainty in the position of the reactions. However, 
this broadening is not sufficiently strong to disrupt the 7?.o profile if the inequality p 1 ^ 2 <C 1 
is satisfied, i.e. if the width of the reaction zone is much greater than the broadening. 
Rather, the fluctuations act to smooth the discontinuous nature of the equally-spaced re- 
actant currents, for which a coarse-graining was necessary in section III. Therefore, from 
equation ( JKj| ) the steady-state intrinsic reaction rate is 

TZi(x r ) = 2cg 2 P{2gx r ) (24) 

MXr) = f 2 f . +0(p 2 ) 

2 cosh {pgxr) 

It also follows that the particle streams have the coarse-grained density distribution seen in 
the ordered case. Hence, figure 2 also represents the profiles seen in the case with Gaussian 
fluctuations in the particle currents. 

The width of the intrinsic reaction zone is ~ (pf?) -1 and from fl23|) in a time AT it will 
typically move (as a whole) a distance ~ (cAT / g) 1 / 2 . Hence, any measurement made of the 
production rate must have AT -C {p 2 cg)~ l . Such a restriction can indeed be satisfied, and 
the production rate that fl2"4"[) predicts can be clearly seen in a Monte Carlo simulation of 
a single system, i.e. the system is self-averaging. 

In the above analysis the steady state was assumed for the intrinsic reaction rate cal- 
culation. However, the arguments also follow through if the time-dependent intrinsic rate 
(|TD|) is used, for p 1. The width of the intrinsic reaction zone Wj therefore varies as 

Wi oc (ct) for t < £/ 
oc (pg)' 1 for t > t T 
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where tj = {peg)' 1 is the relaxation time for the intrinsic rate given in the previous section. 
Hence, at early times the width increases linearly, until finally saturating at a finite time- 
independent width. 

It is unlikely that the present model studied is in the same universality class as the case 
of diffusive reactants. However, it is interesting to note that the asymptotics of the intrinsic 
reaction rates derived for the diffusive case and in the present case of ballistic reactants are 
both of the same form ~ e~ a ^, in the steady state. 

The extrinsic reaction rate 

The extrinsic reaction rate, defined as the reaction rate at (x, t) averaged over all possible 
evolutions, implies an average over all the appropriately weighted initial particle positions 
{y} and {z}. It is therefore equivalent to the sum of probability densities for any (to, n) 
pair to annihilate at (x,t). After the intrinsic reaction rate has relaxed to its steady-state 
limit, i.e. for times t ^> ti, the extrinsic reaction rate can be expressed purely as a function 
of s — (n — to). 

Tl E (x,t) = J2J2P(m,n)G mn (x,t) 

m n 

]imK E (x,t) = VP(s)£ s (:M). 

t— »oo — 

s 

The function Q s (x,t) is the probability density that an (to, to + s) pair meet at (x,t) for 
any value of to with s fixed, in the limit t 3> (cg)^ 1 

oo 2 / j \ 

Q s (x, t) = lim V — exp ( — ((2gx - s) 2 + (2cgt - 2m) 2 ) ) 
*^°° r ~ L 7rTO \4m v ~ ') 

eg 2 (-(2gx-s) 2 \ 

Therefore, the extrinsic reaction rate can be written as a convolution of the intrinsic rate 
over all the allowed paths the intrinsic reaction zone can take (each being weighted by a 
Gaussian). 

TZ E (x,t) 



1Z E (x,t) 

A second characteristic time t E — (p 2 c^)" 1 is now introduced into the system by the Gaus- 
sian noise. This is the timescale on which the uncertainty in the position of the intrinsic 
reaction zone becomes equal to its steady-state width, i.e. the standard deviation of the 
position of the reaction zone's centre is of the order (pg)' 1 . For t t E the extrinsic and 
intrinsic rates are effectively the same, as there is little broadening. However, for times 
t ^> t E the noise in the particle currents dominates, the extrinsic reaction rate becomes 



eg 
y/ncgt 

p 2 cg 3 
V 16vrt 



P(s)ex.p 



■(2gx - 
Acgt 



(25) 



E 



exp 



Acgt 



coslr (ps/2) 



+ 0(p 2 ). 
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a Gaussian and the result || is recovered. This is not surprising as for these times the 
structure of the intrinsic zone is unimportant and can be considered a point-like random 
walker. The width We of the extrinsic reaction rate as a function of time is therefore 

W E oc (ct) for t < ti 
W E oc (pqY 1 for £j < t < £ B 
oc (ct/^) 1/2 for t > t s 

where ij = (pc^) -1 and Ie — (p 2 c£>) _1 . The asymptotic behaviour (large x values) of the 
extrinsic reaction rate also has distinct forms as a function of time 

K E ~ 0(ct - g\x\) for t < tj 
K E ~ e" a|x| for t 7 < t < t B 
7^ ~ e- fel2/t for t > £ s . 

Hence, there is a cross over from exponential to Gaussian behaviour at late times. 



V. DISCUSSION 

In the present work, a new model for the reaction system A + B — > with ballistic reac- 
tants has been introduced and solved exactly. The model includes two types of noise, in 
the reactions (due to the reaction probability being less than one) and in the currents of 
particles (due to disorder in the particles' initial positions). These fluctuations allow some of 
the characteristics of real systems to be exhibited, including a fluctuating A and B overlap 
region where the reactions occur, and a reaction zone that moves stochastically throughout 
the system. Comparison was made between the present model of ballistically moving re- 



actants and that of isotropically diffusing reactants [0 studied under the RG framework. 
In particular, it was noted that the intrinsic reaction rates (the dynamic reaction region 
formed between the fluctuating A and B domains) have the same asymptotic form, though 
it is unclear if the two systems share the same universality class. An important physical 
difference between these two cases is that, in the ballistic model, particle order is preserved 
and interactions between a given pair of particles can occur only once. It would therefore 
be interesting to study the case of reactants that perform biased diffusion, because such a 
system would interpolate between the cases of isotropic diffusion and ballistic motion of re- 
actants. As shown in the appendix, the method outlined in section II can be translated into 
the second-quantisation formalism with the vacancy dynamics described by the SU q (2) al- 
gebra. Using the tools available in this formalism it may be possible to introduce interaction 
terms in the evolution operator that break the order of the particles, thereby introducing a 
bias-diffusive component into the otherwise deterministic motion. 

Another case for further study would be to examine more closely the relation between 
the present model and the partially-asymmetric exclusion process (PASEP) with reflecting 
boundaries - also described by SU q {2) [p4fl . The present model displays short-range cor- 
relations, as does the PASEP with reflecting boundaries. However, the PASEP with open 



boundaries is known to have three phases [25[]-[2J|, one of which has long-range power-law 
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behaviour. It would therefore be interesting to see if the open system also translates into a 
system of reacting particles. 
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APPENDIX. "Second-quantisation" formulation 

As stated in section II, the method used to calculate P(m, n) can be viewed in terms of 
quantum-mechanical interacting particles. Here this formulation is briefly reviewed and the 
recipe for calculating quantities of interest is described. 

The system studied in section II consists of a one-dimensional lattice with sites numbered 
1,2- ■■, etc. A lattice site k, described by the binary variable rjj., can either be occupied 
by a particle or a vacancy. If a vacancy is present at site k then = 1, otherwise rjk = 0. 
Hence, the state of a single system is described by the set of variables 77 = {rjk} and can be 
represented vector in a Fock space, i.e. 



W\v) 



\VUV2, 



Vk---) 



k=l 



Initially the system is full of B particles, represented by the state |0). The evolution of the 
system involves an A particle being 'shot' through the B array at time steps m = 1, 2 ■ ■ -, 
etc. Each A particle passes through a B with probability q, or annihilates with the B with 
probability p. If no B is present at a site the A particle passes through to the next site with 
probability one, see figure 1. After each time step m, the system will have some probability 
Pm{v) °f being in the state rj. The vacancy creation and annihilation operators are now 
introduced 

fc-i 



c+|-..o fc ..-> 



Cfc| • ■ ■ lfc • ■ ■) 



. 1=1 
■ fc— 1 

. i=l 



■■h 



•} 



with the auxiliary relations C£\ ■ ■ ■ lk ■ ■ •) = C&| • • • • • •) = 0. Defining the g-commutator 

as [A, B] q = AB — qBA, the following commutations relations hold for rii < n 2 

(26) 



ri+ (-1+ 



C C 







Cni Cn 



1 



where the final anticommutator is for same-site operators. These commutation relations 
provide a representation of SU q (2), i.e. 'g-deformed' spin-half particles p8| . The evolu- 
tion equation for the system can now be written in terms of these vacancy creation and 
annihilation operators 

\P m ) = Y, P m(v)\v)=f m \0) 
V 



T 



k=l 



where f is the evolution operator. Introducing the left state (S\ = ^'(Vl an Y observable 
A(P m ) represented by the operator A can be written in the following way 

(s\Af m \o) = EJ2(v'\v)A(v)Pm(v) 

(S\Af m \0) = ^A( V )P m (r,). 
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In evaluating such observables it is useful to use a factorisation characteristic of the creation 
operators. If a string of creation operators is arranged, such that each operator acts in 
descending order (with respect to its subscript) on the right empty state, then the following 
property can be used 

(S\C+C+---Cl i \0) = (S\C+ 1 \0)(S\C+\0)---(S\C+jO) (27) 

where the single-vacancy expectations are simply (S\C£\0) = pq k ~ l . 

In the formalism described above, all quantities expressible in terms of operators can be 
evaluated. Therefore 

*K...,n m ) = (S\ ]T C+ 1 C+---C+ n \0) 

perms 

V m (n) = (S\C+C n T m \0) 
P(m,n) = (S\C+f™- l \Q). 

are the forms for each of the objects calculated in section II. These can be evaluated by using 
the commutation relations ( p6|) in the following way. Firstly for a given string, commute 
all annihilation operators to the right of the string. These objects when acting on the left 
state (5*| will leave strings containing only creation operators. These strings can then be 
rearranged by using the appropriate relation in (EB[) so that they act on the zero state 
in descending order of subscript. Finally the factorisation property (|27j) can be used to 
calculate the expectation value. 

Interestingly, the algebraic structure outlined here is identical to that used in the de- 
scription of the partially-asymmetric exclusion process (PASEP) with reflecting boundary 
conditions, therefore corresponding to a noisy Burgers equation with zero average current 



24]. Microscopically, the PASEP describes systems of particles that hop in a preferred 
direction with a repulsive interaction. Assuming that the particles hop with a leftwards 
bias, the density at site n seen when m such particles are confined to a semi-infinite lattice 
([1 : oo]), is identical to the vacancy density V m (n) given by equation (|T0|). 
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Figure captions 



FIGURE 1: 

(i) A realisation of the model described in section I. The solid lines are paths of particles 
in space and time, with the dotted lines the trajectories annihilated particles would have 
taken. The A and B particles start at positions y m and z n respectively and move with 
fixed velocities va = c and vb = —c. An intersecting pair of particles either annihilates 
(probability p) or continues unaffected (probability q = 1 — p) . 

(ii) The simplified, mapped system described in section II. At unit time steps A particles 
are 'shot' through the B lattice, passing through each occupied site (with the same micro- 
scopic reaction probabilities as above). A reaction occurring at a site produces a vacancy, 
and it is the statistics of these vacancies that is used to derive P(m,n), the pair-reaction 
probability. 

FIGURE 2: 

(i) The steady-state density profiles of the A and B reactants at x r = s/g for the case 
p = q = 1/2. The densities given in equation (|17D involve the vacancy densities (|10|) for the 
exact case (circles and squares) and ([□]) for the 0(p) approximation (dashed lines). 

(ii) The steady-state reaction rate at position x r = s/g, again for p = q = 1/2. The 
the exact result (circles) is given in equation (|T6|) with x r measured from the origin, or in 
equation ( p4|) with x r measured from the centre of the fluctuating overlap region. Also plot- 
ted are the mean-field (triangles) and the 0(p) approximations (dashes) given in equations 
(p0|) and (|19| ) respectively. Comparison between the exact and mean-field results show the 
reaction rate to be most correlated at the centre of the overlap region. 
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(i) Particle World Lines. 
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(i) Reactant Density in the Steady State 
for the Intrinsic Reaction Zone. 




-10 -5 5 10 



distance from overlap centre 



(ii) The Intrinsic Reaction Zone 

in the Steady State. 
— i 1 1 1 1 1 1 — 

O — O^/ reaction rate 
A — fonean field 




distance from overlap centre 



Figure 2: 



20 



